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We build the first analytic empirical potential for the most deeply bound Li 2 state: b . Our 

potential is based on experimental energy transitions covering u = 0 — 34, and very high precision 
theoretical long-range constants. It provides high accuracy predictions up to v = 100 which pave 
the way for high-precision long-range measurements, and hopefully an eventual resolution of the 
age old discrepancy between experiment and theory for the Li(2^5) -I- Li (2^P^ Cs value. State of 
the art ab initio calculations predict vibrational energy spacings that are all in at most 0.8 cm~^ 
disagreement with the empirical potential. 

PACS numbers: 02.60.Ed , 31.50.Bc , 82.80.-d , 31.15.ac, 33.20.-t, , 82.90.-|-j, 97, , 98.38.-j , 95.30.Ky 


There is currently a rather large discrepancy between 
the best atomic 3e“ ab initio calculation [l|, and the 
most current empirical value ii, for the leading long- 
range Li(2^S') — Li(2^P) interaction constant (Ca), de¬ 
spite the latter being the most precise experimentally 
determined oscillator strength for any system, by an or¬ 
der of magnitude [J. Li 2 (&, l^IIii) is one of the molecular 
states that dissociates to Li(2^5')-|-Li(2^P), and therefore 
its long-range potential has this interaction constant. 
This “&-state” is also the deepest Li 2 potential, and out 
of the five lowest Li 2 states, the 6-state is the only one 
for which an analytic empirical potential has never been 
made. 

Since the highest 2000 cm“^ worth of vibrational levels 
of the 6-state still have not been observed, and part of 
this region is now accessible by current ultra-high preci¬ 
sion PA (photoassociation) technology Bi, an analytic 
potential would be very useful for making predictions to 
assist in observing the missing levels. The 6-state also 
mixes strongly with the A(l^Eii)-state , which has by far 
the most precisely determined excited molecular state po¬ 
tential in all of chemistry, yet still has a rather large gap 
of missing data in the middle of its energy range [3, Q . 

Finally, the 6(l^n„)-state has been a key doorway to 
the triplet manifold, and was directly involved in the 
measurements for a vast number of other triplet states 
such as 23Ag [|, S^Hg 0, Bl, 23E+ B [a, 

2^ng Bfllill: l^S-, and other un¬ 

determined 3A states [H, [l^ . Some of these more highly 
excited triplet states (namely 33 E+, 23 E+, 23ng, and 
l3Ag) are so thoroughly covered by these spectroscopic 
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measurements, that global empirical potentials can be 
built for them too. For this, an analytic potential for the 
6-state would be used as a base. 

In this work we will build analytic empirical potentials 
for the 6-states of all stable homonuclear isotopologues 
of Li 2 . Previous work has shown that analytic empirical 
potentials for the c(l3Eg)-state were able to predict en¬ 
ergies correctly to about 1 cm“^, in the middle of a gap 
of > 5000 cm“3 where data were unavailable @,[13 , and 
this was much better agreement than was obtained with 
the most sophisticated Li 2 ab initio calculations of the 
time [T^ . 

It was recently shown that the best ground-state rota¬ 
tionless ab initio potentials for the 5e“ molecules BeH, 
BeD, and BeT, were able to predict vibrational energy 
spacings to within I cm“^ for all measured energy levels 
except one. The 6-state of Li 2 might be expected to be 
more challenging ab initio because it (I) is an excited 
state, (2) has one more e“, and (3) involves many more 
vibrational energies. We will therefore compare our an¬ 
alytic empirical potentials for the 6-state of 3’®Li2 and 
^dLi 2 with the most state-of-the-art ah initio calcula¬ 
tions, which were published recently in [l^ . 

Table U summarizes all experiments we could find 
which provided information on rovibrational levels of the 
6-state . Unfortunately attempts to recover the data from 
B in, [20l - [^ were unsuccessful, but we were still able to 
include all data from the other experiments in our study. 
Furthermore, it is noted that the 6-state was also involved 
in various other studies but these just made 

use of rovibrational levels that were already determined 
in the studies listed in Table [H in order to access levels 
of other electronic states. 
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Table I. Summary of experiments involving Li2(b, l^n^) 


Isotopes Year 

Type 

States Involved 

Unc. (cm ^) 

V 

^ Data Included in dataset 

Source 


1985 

LIT 

A(HE:),b(i3n„) 

a 

8-25 

100s 


Preuss &: Baumgartner 


1985 

OODR 

(MS^Ag) ^ b(l3n„) 

1 

0 

100s 

- 

Rai et al [T^ 


1992 

PFOODR 

A(iiE+),b(i3no 

? 

19 

3 

- 

Li et al. l26ll 


1996 

cw PFOODR 

eAg) ^ Ki^no 

? 

0,1,4 

7 

- 

Li et al.|l6| 


1996 

GIF 

(i^A,,) ^ b(i3n„) 

? 

0-11 

7 

- 

Weyh et al.|22l| 


1997 

PFOODR 

(23n9) ^ fe(i3n„) 

0.005 

1-27 

178 

/ 

Russier et al fTl| 


1997 

PFOODR 

(i^A,,) ^ b(i3n„) 

0.005 

1-25 

234 

/ 

Russier et al fTl(l 


1997 

GIF 

(i^Ag) ^ b(i3n„) 

0.005 

0-7 

314 

/ 

Russier et al fTl| 


2001 

cw PFOODR 

A(FEi),b(i3n„) 

? 

15,22 

2 

- 

Lazarov Sz Lyyra [Toll 


1985 

LIF 

A(FE+),b(i3n.) 

a 

8-20 

100s 


Preuss &: Baumgartner [25II 

®'®Li2 

1983 

GIF 

(S^Ag) ^ 6(l3n„) 

0.3 

0-10 

240 

- 

Engelke &: Hage 


1985 

LIF 

A{PE+),b{i^n^) 

a 

8-18 

100s 


Preuss &: Baumgartner |25|I 


1985 


A(iiE+),b(i3n„) 

7 

9,15 

2 

- 

Xie & Field 113 , 111 ] 


1986 

PFOODR 

(23n3) ^ b(i3n,,) 

0.5 

0-17 

-170 

32 lines recovered 

Xie & Field 0 


1986 

PFOODR 

(i^A,) ^ b(i3n„) 

0.02-0.13 

0-11 

7 


Rice, Xie & Field [23,[30ll 


1988 

GIF 

(i^A^) ^ b(i3n„) 

0.2 

0-11 

-200 

- 

Schmidt et al. \2iJ 


1992 

GIF 

(l^Ag) b(l3n„) 

0.003 - 0.07 

0-9 

599 

/ 

Linton et al. [l2ll 


TOTAL 


0.003 - 1 

0-27 

7 

1357 

mini 


“The measurements were on a-levels of the T-state, and information about the b-state n-levels that perturbed those T-state levels was 
inferred indirectly. 


I. HAMILTONIAN 

The rovibrational energy levels and wavefunctions for 
isotopologue a with reduced mass are treated as 

I 

(1 + 5«(r))^ V'v.wW = K.iv(r)V'j;.iV • (1) 


the eigenvalues and eigenfunctions in the effective radial 
Schroedinger equation: 


Here Va{r) and gair) represent the “adiabatic” poten¬ 
tial and the “non-adiabatic” rotational 5 -factor. The 
adiabatic potential can be represented as a “Born- 
Oppenheimer” potential (which is mass-independent), 
plus a (mass-dependent) shift due to the diagonal cor¬ 
rection to the Born-Oppenheimer approximation: 


V^{r)=VBo{r)+AV^{r). ( 2 ) 


The AVa(r) correction can be approximated by the ex¬ 
pectation value of the nuclear kinetic energy operator in 
the molecular electronic wavefunction basis 
For homonuclear diatomics it is given by 



Qair) = -7 


Po^ir) = 


2g,c 


1C,a) + A2Va(r) 

(3) 

r) + P^iR) + S^iR) 

(4) 

/ 52 \ 



(5) 

- {^l+^Da 

(6) 
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-|;(i:v.v,) ( 10 ) 


= ^^Bo(?’) + r-^VBo{r)^ + AS'a(rX.ll) 

where z represents the internuclear axis, and Ly are 
then projections of the total electronic orbital angular 
momentum, i,j represent indices for individual electrons 
of the molecule, and the first term of Sair) has been ex¬ 
pressed in terms of the average electronic kinetic energy 
{Te,a) and then re-expressed in terms of k^o(?') using the 
virial theorem . We can define a long-range term 

^LR,a(?') by evaluating Pa{r) in the long-range Heitler- 
London basis, where electron overlap is zero. Pn(r) is 
then expressed as Plr.o (?’) plus a correction AP„ (r) : 

I 


A14(r) 




^{yBo(r)+rlvso(r) 




me 

4/ra 


VBoir) + r—VBo{r) 


where me is the electron mass and Mk is the mass of 
the constituent nucleus of the molecule. Note that 
until now, the terms containing Sk{r) represented the 
entirety of Eq. [T3 so less of AVa (r) was described by 
theoretically known expressions, and more was 
described by empirical fitting functions [s^. 

The only part of the Hamiltonian in Eq. [T] that is 
missing is now the non-adiabatic term ga{f)- This is of¬ 
ten represented by model functions Rk(r) for each atom: 

ga{r) . (18) 

k 


II. EMPIRICAL POTENTIAL AND 
BORN-OPPENHEIMER BREAKDOWN (BOB) 
CORRECTIONS 


We now wish to determine empirical functions for 
Vbo(?'), S{r), and R{r) that accurately reproduce all 


Pa(r) =PLR.a(r) + AP„(r) (12) 

PLR,a(r) = ;^^Zfe(4 + l), (13) 

z/rQ,r ^ 

where Ik represents the orbital angular momentum of the 
electrons in constituent atom k of the molecule. Herein 
we restrict our attention to the Pstate of Li 2 which dis¬ 
sociates into Li(5') -I- Li(P): 


PLB.,a{l") — -2 + 1) + Ipilp + 1)) (14) 

ZllciT 

2h^ 

= -2 ’ =0,lp= 1) . (15) 

While we know that in the long-range limit, Qair) will 
become a constant [3l|, APa{r) will be zero, and ASa{r) 
will be small , no other information about these terms 
is known. Therefore, we may re-write the diagonal Born- 
Oppenheimer correction (DBOC) in terms of what we 
know, and then represent these parts that we don’t know, 
by model functions Sk(r) for each constituent atom k of 
the molecule: 


+ (A5„(r) + AP„(r) + Qa(r) + A2Vc{r)) (16) 


measured energies when using the Hamiltonian of Eq. 

m 

There is a gap of more than 2000 cm“^ (> 60 THz) in 
experimental information between the highest observed 
level of Li 2 (6, l^Hii), and its dissociation energy. This 
means that when building an empirical potential that 
aims to be relevant in the large data gap, it is very im¬ 
portant to take great care in ensuring the potential be¬ 
haves physically correctly in the extrapolation region. In 
2011 the MLR (Morse/long-range) model was used in a 
fit to build empirical potentials from spectroscopic data 
for the c (l^E+) —states of ®’®Li 2 and ^’^Li 2 , where there 
was a gap of more than 5000 cm“^ between data near the 
bottom of the potential, and data at the very top 0. In 
2013 spectroscopic measurements were made in the very 
middle of this gap and it was found that the vibra¬ 
tional energies predicted by the MLR potential from 0 
were correct to within about 1 cm The present case 
for the 6-state is in some sense more interesting because 
there is no data at the top helping to anchor the potential 
with the right shape near dissociation. 
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Figure 1. When considering the data for all isotopologues, levels up to n = 32 have been observed. For ®’®Li2, levels n = 33 — 40 
are accessible W Kirk Madison’s group in University of British Columbia with the laboratory and method used in studies of 
the A d and c (l^E+) [| states. The remainder of the levels have not been observed. 



However, as for the c-state, the MLR model is still 
expected to be able to represent the physics in the ex¬ 
trapolation region faithfully since the correct theoretical 
long-range is built into the model. Having this long-range 
physics accurately built into the model is almost as help¬ 
ful as having data in the long-range region, as was the 
case of the c-state. MLR-type empirical potentials have 
now successfully described spectroscopic data for many 
diatomic [2, [1, [2, [13, and polyatomic [H, 

systems. Therefore, we will proceed to use the MLR 
model to describe Ubo(’')- 
The MLR model is defined by 


, (19) 

V w(re) J 

where De is the dissociation energy, Cg is the equilibrium 
internuclear distance, and the polynomial /3(r) is 

N/3 

/3(r) = + (l “ ^ Pi , 


/3oo = lim /3(r) = In 

r—^oo 



( 21 ) 


Equations [H] and [201 also depend on the radial variable 


yn’’'{'r) = 




ref 


( 22 ) 


' ref 


where the reference distance Cpef is simply the equilibrium 
distance in most cases, but can be adjusted to optimize 
the fit to equation [H 

It is well known [2, [ij] that for large r we have 


U(r) ~S)e-M(r) + --- , (23) 

therefore the long-range behavior of the potential is de¬ 
fined by u(r), and the short to mid-range behavior is de¬ 
fined by /3(r). In the b (l^n„) state, a spin-orbit interac¬ 
tion emerges at large internuclear distances, which splits 
the potential into four components. Therefore, four dif¬ 
ferent potentials (r) can be defined to have the same 
/3(r) defining the short-range behavior where there is no 
significant splitting, and to have four different u 
defining the long-range where the splitting occurs. 

For large r where the spin-orbit interaction becomes 
strong, it is dangerous to label the spin angular momen¬ 
tum S and orbital angular momentum A separately, as 
in the molecular term symbol Instead, these two 

momenta are combined into a total electronic angular 
momentum H. For S' = 1 and A = 1, so there are 
states with four possible symmetries in the H representa¬ 
tion: 0+, 0~, lu, and 2„. Each of these four states has a 
slightly different behavior at large internuclear distances. 
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due to coupling with states that have the same symme¬ 
try in the representation, but different symmetry in 
the A representation. This coupling has been described 
in [ 6 l| and has been used for building appropriate an¬ 
alytic empirical potentials for the A(l^S+) [ 3,0 and 
c(liS+) I, [13 states of Li 2 . The long-range function 
u{r) is defined separately for each spin-orbit state: 



(r) 

, 0+ 


(r) 

,0- 


(r) 

; lu 

y 6 . 2 „ 

(r) 

2 

5 


(24) 


Each of these functions is an eigenvalue of a matrix 

I 


for each state. These matrices are given in the subsec¬ 
tions below, in terms of the (positive) spin-orbit splitting 
energy AE, and neglecting exchange interaction terms. 


A. The Oj states 


In addition to the 6 (l^n„) state, the other A state that 
can give rise to 0 + symmetry is the A( 1 ^E+) state [^ . 
The interstate coupling is therefore given by the 2x2 
matrix [ 6 l[: 


(r) = 


/ a(i1s+) 6(i^nu) 

1 '+2Cr^ ^ 

3 

m=3,6,8 

^ V ’-C^ > 

3 r™ 


# E 


c. 




m=3,6,8 

9,10,11,... 


a(i1s+) b(i3n„) 

Aif+i E 


V 


m=3,6,8 

9,10,11,... 


m—3,6,8 
9,10,11,... 


(25) 




where the lower energy eigenvalue (r) comes from A (l^E+) and approaches the dissociation limit of Li( 2 ^S'i/ 2 ) + 
Li ( 2 ^Pi/ 2 ), and the higher energy eigenvalue (r) comes from b (l^n„) and approaches the dissociation limit of 
Li(22S'i/2) -hLi(22p3/2) [el. 


B. The Ou states 

In addition to the 5(l^n„) state, the other A state that caii give rise to 0~ symmetry is the 2a(2^E+) state [63 |. 
The interstate coupling is therefore given by the 2x2 matrix [63 : 


u°~(r) = 


/ 2a(2^S+) l>(l='n„) 

/ 1 C„A '+2CA ’ 

3 

m=3,6,8 
9,10,11,.. 

# E 

m—3,6,8 
\ 9,10,11,... 


C, 


2.(23^+)_ 


___ 2a(23E+) 6(l3n„) \ 

Z 1 '-cA ’ \ 

3 Z^ r3 

m=3,6,8 
9,10,11,... 

-AE+I ^ 1 


m—3,6,8 
9,10,11,... 


(26) 




where the lower energy eigenvalue comes from b (l^n„) and approaches the dissociation limit of Li( 2 ^S'i/ 2 )-l-Li {2^Pi /,^, 
and the higher energy eigenvalue comes from 2 a ( 2 ^E+) and approaches the dissociation limit of Li( 2 ^S'i/ 2 )-l-Li ( 2 ^^ 372 ) 

M- 


C. The lu states 


In addition to the 6 (l^nn) state, the other A states that can give rise to I„ symmetry are the 2a (2^S+) state and 
the B (l^n„)state [63|. The interstate coupling is therefore given by the 3x3 matrix [^ : 


ui“(r) = 


(27) 
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(28) 


where the lowest energy eigenvalue comes from b (l^n„) and approaches the dissociation limit of Li( 2 ^S'Y 2 ) + 
Li( 2 ^Pi/ 2 ), and the middle and highest energy eigenvalues come from S (l^n„) and 2 a ( 2 ^E+) respectively and 
both approach the dissociation limit of Li( 2 ^S'Y 2 ) + Li |^ . 


D. The 2u state 


The 2„ state approaching the dissociation limit of 2S + 2P is alone in its symmetry, and approaches Li( 2 ^S'Y 2 ) + 
Li ( 2 ^P 3 / 2 ) It has the long-range function 


‘W = - al;- ^ 


a 


b(i^n„) 


\ 


m=3,6,8, 

9,10,llv 


/ 


= - \ AE- 


a 


b(i"n„) ^b(i^n„) ^b(i3n„) ^b{Pa^) 


C, 


Co 


b(i3n„) b(i^n„) 

Oio 


.10 


nil 


(29) 


(30) 


E. All four flu states combined 


F. Quadratic corrections and damping functions 


One can imagine an experiment which obtains spec¬ 
troscopic measurements for all of the four states, and 
fits to all of this data simultaneously by using the appro¬ 
priate eigenvalues of the 8 x 8 matrix below: 




\ 


(31) 


oi 

oj 

11 

'*^12 


ot 

'21 

U22 


0,7 

07 

^11 

^12 

07 

07 

U21 

'^22 


V 


\ 


MlJ 

U 12 

'^13 


lu 


■“21 

^22 

■^23 

'^31 

'^32 

"'^33 


.(32) 


However, Fig. [T] shows the data region, and Fig. [2] 
shows that the spin-orbit splitting does not seem to be¬ 
come apparent until well past this region. Since the mea¬ 
surements that have been done on the &-state thus far 
are far away from the effect of the spin-orbit splitting, 
we choose to use the simplest spin-orbit long-range func¬ 
tion: u^’^"(r). 


/ >,2 

Since the leading term not shown in Eq. [23] is , 
the contribution of the C 3 terms to the long-range form 
of the potential, will interfere with the desired Cq and Cs 
terms, and all Cg and Cn terms will therefore have spu¬ 
rious contributions from the cross-terms formed by the 
products of the C 3 terms with the Cg and Cs terms re¬ 
spectively. We fix this in the same way as was done for Cq 
and Cg in iii, [TtI . , by applying a transformation 
to all Cg, Cg, and this time also Cn terms: 

Cg ->■ Cg -I- (33) 

Cg ^ Cg + (34) 

Cii ^ Cn + (35) 

where the transformation in Eq. [33] has to be made 
first due to Eq. [SUs dependence on Cg. 

Additionally, The long-range formulas in terms of Cm 
constants in the above sub-sections were derived under 
the assumption that two free atoms are interacting with 
each other, and there is no overlap of the electrons’ wave- 
functions as in a bound molecule. To take into account 
the effect of electron overlap, we use the damping func¬ 
tion form from 
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Figure 2. In the long-range region, the splitting of the b state into four spin-orbit components becomes increasingly 

obvious. This is well beyond the region where data is available, so our empirical potential fits to a model with the simplest 
long-range potential energy function (the 2„ state, since there are no other 2^ states that are nearby in energy, see Eqs. I29I32II . 
The fine-structure splitting of 0.33532461313 cm“^ comes from measurements in [^ . 



Cm ^ CmD\:} (r) (36) 



where for interacting atoms A and B, p = pab = 
in which px = is defined in terms of 

the ionization potentials of atom X, denoted (7^), and 
hydrogen(/^). We use s = — 1, which as shown in [i^ . 
means that the MLR potential in Eq. [12] has the phys¬ 
ically desired behavior V oc i/r^ in the limit as r —>■ 0. 
For s = — 1, the system independent parameters take the 
values = 3.30, and = 0.423 [d^. 

G. Long-range constants 

In previous studies of the A(1^I1+) state d, and 
c(13E+) state @,[13 , it was found that the most precise 
theoretical values of C 3 known at those times [^, did 
not fit as well with the measurements of the high-lying 
vibrational levels near the dissociation, as the values of 
C 3 obtained by setting it as a free parameter determined 
by a least-squares fit to the data. 

However, for the 6 (l^n„)-state, no measurements of 
such high-lying vibrational levels have been made, so 
such an “empirical fit” to C 3 is impossible, and we will 
have to use the most precise theoretical value known. 


For ^7Li2 this is the value from and for ®’®Li 2 this is 
an unpublished value from Tang calculated in 2015 [^ . 
These values are listed in Table [iTl along with the the¬ 
oretical values for the higher-order Cm constants used 
in our analysis (it has not yet been possible to fit these 
higher-order to > 3 constants to spectroscopic data in 
any direct-potential-fit analysis, so they are held fixed). 
For TO > 8, no finite-mass corrections have been calcu¬ 
lated yet. 

H. Dissociation energy De 

At the time of carrying out our analysis, the best ex¬ 
perimental value for Tie of which we were aware, was the 
1983 value from 12145±200 cm“^. In a recent study 
on BeH [53, the gap between the highest observed level 
and the dissociation asymptote was ^ 1000 cm“^, and 
the fitted value of Dg varied by about 400 cm“^ as pa¬ 
rameters such as Tref, p, and q were changed. For the 
present case of the 6-state of Li 2 , the data region stops 
more than 4000 cm“^ below the dissociation asymptote, 
so we do not expect to be able to determine Se^'iiy more 
precisely than the 1983 experimental value. However, we 
still tried, by letting Te be a free parameter, and we in¬ 
deed found that the fitted values varied by more than 
400 cm“^. Therefore, it might make sense to use exper¬ 
imental value from [20|| which was claimed to be within 
400 cm“^. 

However, it is expected that the ab initio value from 
[13 correct to within much less than 400 cm“^. This is 



















Table II. The best currently available long-range constants and their sources (in Hartree atomic units). °°Li2 represents a Li2 
molecule where both nuclei have infinite mass, since finite mass corrections have not yet been calculated for Cm coefficients 
with m > 9. 



A(1 

®'®Li2 

Is,) 

^’^Li2 

2a (23E„) 

®'®Li2 ^’"Li2 

B (Pn„) 

®’«Li2 ^’^Li2 

b(i"n,) 

«’®Li2 ^’^Li2 

Ref. 

C3 

11.0009 [H 

11.0007 [Hi 

-11.0009 [H -11.0007 IH 

-5.5005 IH -5.5004 H 

5.5005 H 5.5004 |H| 

- 

Ce 

2076.19(7) 

2076.08(7) 

2076.19(7) 2076.08(7) 

1407.20(2) 1407.15(5) 

1407.20(2) 1407.15(5) 


Cg 

274137(6) 

274128(5) 

991104(5) 991075(6) 

48566.9(4) 48566.4(2) 

103053(1) 103052(1) 

m 


OO 

Li2 

°°Li2 

“Li2 

°°Li2 


Cg 

2.2880(2) X10® 

-2.2880(2) X 10® 

-5.173(1)x104 

5.173(l)xl04 

i 

Cio 

3.0096x10'^ 

-1.2113x10® 

8.9295x10® 

9.1839x10® 

0 

Cii 

-5.930 X 10'^ 

-5.0321x10® 

-9.924 X 10’^ 

2.652x10'^ 

\i 


Figure 3. Various (A(a)p(q^ models with the same quality of fit are in agreement with each other, until large values of r at which 
the Np = 7 model dips significantly below the red curve representing the theoretical long-range behavior. This dip was obvious 
in the —r^V(r) vs plots for all models mentioned in the text except for the three Np = 11 cases shown. The inset shows 
—r^V{r) for the four mentioned models, subtracted from the red theoretical curve. All models match the theory for very large 
r (small but the cases with higher p and/or q match the theory for a much larger range of r. Each green circle represents 

an outer classical turning point of the ab initio potential. In the inset, the ab initio values lie too high above 0 to be seen. 



l/r3 [A-3] 


because we systematically checked the ab initio 'De values 
for all electronic Li 2 states calculated in , and found 
that they were at most 68 cm“^ different from the best 
experimental value, even when the experimental values 
were known to as high of a precision as 0.0023 cm“^ (see 
Table Hill) . Furthermore, the ab initio value for the 6-state 
was within the 400 cm“^ conhdence interval given by the 
1983 experimental value [2^ discussed in the previous 
paragraph. Therefore, we decided to fix our 2)^ value at 
the ab initio value of 12166 cm“^ and to only allow the 
other parameters be free parameters for the remainder of 


the Htting analysis. 


After the completion of this work, we discovered that 
a much less known paper co-authored by one of the same 
authors from [2^, reported a more precise De value of 
(12180.6±0.6) cm“^ just over 4 years afterwards 2j|, but 
it is not clear in the paper how this value, nor its uncer¬ 
tainty is obtained. Particularly, it is not clear whether 
this is a purely empirical value, or if it also uses the ab 
inito potential which is part of the subject of the paper. 
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Table III. Dissociation energies 0e in cm“^ of various electronic states of Li 2 . “obs - calc” denotes the ab initio value subtracted 
from the empirical value. Numbers in parenthesis indicate the experimental uncertainty in the last digit(s) shown. Where 
possible, ®’®Li 2 values were used, for the empirical values. When empirical values for the Se of ®’®Li 2 were not available, the 
empirical value for ^’^Li 2 is given. All ab initio values are for °°Li 2 which is a hypothetical Li 2 molecule with infinite mass. 


2S + 2S 

X(UE+) a(UE+) 





ab initio [19] 

8466. 334. 





empirical 

8516.7800(23) [2] 333.7795(62) [5] 





obs - calc 

51. 0. 





2S + 2P 

a(Ue+) R(Un„) 

2 V(2iE+) 

c(Ung) 

6 (i®na) 

c(UE+) 

ab initio [f^ 

9356. 2930. 

3289. 

1426. 

12166. 

7080. 

empirical 

9353.1795(28) [2] 2984.444(110) [71] 3318(66) [72, 73] 1422.5(3) [74] 12145.(200) [20] 

7092.4926(86) [5] 

obs - calc 

-3. 54. 

30. 

-3. 

-21. 

12 . 

2S + ZS 

3V(3^E+) (1st min.) 

2 A(UE+) (1st min.) 

2 A(UE+) (2nd min.) 

ab initio [1^ 

8290. 

5608. 

5389. 

empirical 

8317. [75, 761 

5615. 

[77, 7^ 

5321. 

[77, 7^ 

obs - calc 

27. 


7. 

-f 

38. 

2S + W 

4V(4iE+) (1st min.) 

2B{2^n^) 

2C{2^\1,) 

D(UAg) 

2d{2^llg) 

ab initio [1^ 

8380. 

6481. 

7773. 

9592. 

8505. 

empirical 

8349. [79] 

6455. [^ 

7773.7(3) [781 

9579. [^ 

8484. [^ 

obs - calc 

-31. 

-26. 

1 . 

-13. 

-21. 


I. Choice of model parameters 

Using the Hamiltonian of Eq. (TJ we fit the parameters 
of Eq. [Tni to the 1234 data, with the involved energy 
levels of the upper states 2^ng and l^Ag treated as free 
parameters. All fits to the data were done using the freely 
available program DPotFit 2.0 [s^. Starting parameters 
for the fits to Eq. [11] were found by fitting to an RKR 
potential using the freely available program betaFit 2.1 
1^ . The RKR potential was made using the program 
RKRl 2.0 using the Dunham coefficients found in Ta¬ 
ble IV of [l2||. 

The quality of a fit was determined by the 
dimensionless root-mean-square-deviation (dd) which 
scales each deviation between an energy predicted by 
the model (Ecaic) and the corresponding measurement 
{Eohs), by the uncertainty of the measurement (uobs), 
for all iVdata measurements: 


dd = 


1 -^data 

1 ^ 

'^-^calc(^) -^obs(^) \ 

fVdata 

V Uohs{.i) / 


(38) 


In previous studies of the A(1^E+) state and 


c(UEg) state i , it was determined that there was no 
benefit in including long-range terms beyond Cg, because 
the data for the high-lying rovibrational energies began 
to deviate from the theoretical long-range potential en¬ 


ergy curve at distances shorter than the distance where 
Cg began to give a noticeable effect on the long-range 
function ULR(r) (see Fig. 6 of Q for example) . 

However, for the 6-state, no data exists in the long- 
range region, so it might make more sense to include 
more Cm terms in ULR(r) in order to anchor the poten¬ 
tial somewhat appropriately in the >2500 cm“^ gap at 
the top of the potential well where no data exists to guide 
the potential. Nevertheless, we first followed the A— and 
c—state studies and only used up to Cg. We found an 
excellent fit with dd = 0.95 with only A(a = 7, p = 6, 
q = 2. However, the long-range behavior of this poten¬ 
tial was in vast disagreement with the long-range behav¬ 
ior expected by theory (see Fig. El). This is because with 
(P:?) = (6,2), the long-range form of the potential de¬ 
scribed in Eq. [23] does not “turn on” until too high a 
value of r (a larger value of r is needed for (r) and 
(r) to become sufficiently close to their limiting val¬ 
ues of 1. 

We can often encourage the the long-range form to 
“turn on” earlier by increasing p and/or q, which often 
comes with the expense of requiring a higher polynomial 
degree Np to recover the same dd. We explored models 
with p G {6, 7,8,9} and q G {2,3}, including the Cg 
term in mrr for p > 7, Cio for p > 8 , and Cn for p = 9. 
Ideally we would always use as m any Cm constants as 
are known, but as explained in [^, [l7i | , the value of p in 
Eqs. Eni and [201 need to satisfy p > wiast — Wfirst where 
















10 


TOiast and mgrst represent respectively the last and first 
Cm terms included in ui^E.{r). We also note that C 12 for 
the Li(2^5') +Li (2^P) asymptote is not available, as far 
as we know. 

We found that if p < 9 and/or g < 3, the long-range 
behavior does not turn on until about r = 32 A (in the 
very best cases), while the m-dependent Le Roy radius 
[851 ca lculated from the radial expectation values found 
in ^ suggests that the long-range behavior should turn 
on before r = 10 A. With {p,q) = (9,3), we found a 
fit with A/3 = 11 and dd = 0.94, where the long-range 
behavior turns on at about r = 20 A (see Fig [3)). In¬ 
creasing g to 4 would likely turn the long-range behavior 
on at closer to the m-dependent Le Roy radius, but no 
(Pi?) = (9; 4) fits with Np < 11 had a, dd < 1 and we 
needed to push to Np > 13 in order to match the dd of 
the best {p,q,Np) = (9,3,11) fits. Using such a high- 
degree polynomial, when the data only required Np = 7 
for a good fit, can be dangerous in terms of the poten¬ 
tial’s extrapolation in the regions neither constrained by 
data nor built-in Cm constants. In this respect, we also 
tried (p, q) = (9,3) fits with Np = 10 for various r^ef 
values, but no such fit had a dd < I. 


J. Born-Oppenheimer breakdown (BOB) 
corrections 

With our best MLR model: MLRp™^ {Np) = 
MLRgJ (II), we attempted to add adiabatic {S{r), from 
Eq. [T71) and non-adiabatic {R{r), from Eq. fTSl) BOB cor¬ 
rections with the same model functions as used in pre¬ 
vious studies of LE since these models were improved 
in 2009 [ 3 , i S [s^. It was surprising that despite 
there being 599 ®Li 2 data (with Umax = 9, A^max = 46) 
and 696 ^Li 2 data (with Umax = 27, A^max = 27), adding 
BOB correction functions did not improve the fit. Even 
when fitting to 3 adiabatic BOB parameters and 3 non- 
adiabatic BOB parameters, the dd went down by less 
than 1%. This is unexpected when there is just as much 
data for each isotopologue, and there is such a big dif¬ 
ference in the highest u and N levels observed for each 
isotopologue. 

Nevertheless, it seems that the isotopologue shifts due 
to the kinetic energy term in the Hamiltonian, and due to 
the mass-dependent BOB corrections incorporated from 
Eqs. [3] and [TH are the only significant sources of en¬ 
ergy difference between ®Li 2 and '^Li 2 for the 6 -state 
(within our data’s precis ion) . This may also explain why 
the ab initio potentials [l9| calculated assuming an in¬ 
finite molecular mass managed to predict both the ®Li 2 
and ^Li 2 energies so fabulously (see discussion in Section 
EL] and Table IV), while the ab initio BOB correction 
functions for the not much lighter molecule BeH, were 
so crucial in matching the ro-vibrational energies pre¬ 
dicted from the ab initio and empirical potentials |54l . l87| . 
Therefore, the final potential that we recommend, which 


is the same for both ®Li 2 and ^Li 2 in the 6 -state except 
for the small mass-dependent contributions coming from 
the kinetic energy and the un-colored terms in Eq. 1161 
does not contain any empirically fitted S{r) and R{r) 
BOB correction functions. 


K. Sequential rounding and re-fitting (SRR) 

Observing the predicted values for yielded by 139 
different fits which had dd < 0.957 (within 1.5% of the 
optimal fit, which had dd = 0.942), we see that no fit 
predicted an outside the range (2.589 825 < Vg < 
2.589 871) A, regardless of the values of {p,q,rge{,Np), 
though the more extreme predictions of Vg within this 
range corresponded to fits with Np < 9. Based on this 
observation, we recommend the value rg = (2.589 848 ± 
0.000 023) A, which is the average of these upper and 
lower bounds, with the uncertainty being the distance 
from the average to either bound. 

We then re-fitted the potential to the data, but with 
re fixed at 2.589848 A, once with the DPotFit setting 
I ROUND = —1 and once with I ROUND = -|-1 in order to 
implement the SRR procedure described in [s^ and in 
the DPotFit manual [82] ■ Neither of these cases affected 
the 3-digit value dd = 0.942. The IROUND = —1 fit ended 
up with 2 more total digits than when IROUND = -|-1 was 
used, but had a lower dd in the 4th digit, and has the 
more elegant feature that the number of digits in fdi de¬ 
creases monotonically with increasing i. Therefore, we 
recommend the potential with IROUND = — 1, whose pa¬ 
rameters are listed in Table IlYl 


L. Vibrational energy spacings of the 
recommended Li 2 (b,® Hu) potential, and comparison 
to best ab initio potential 

Very recently, a review paper on the 5e“ systems BeH, 
BeD, and BeT [sj revealed that the state of the art 
ah initio potentials [s^l (which used MR-ACRE/aug-cc- 
pCV7Z(i), a further estimate of electron correlation ef¬ 
fects beyond the approximations of MR-ACPF, second- 
order DKH scalar relativistic corrections, and mass- 
dependent BOB corrections), predicted vibrational en¬ 
ergy spacings with up to at most 1.8 cm“^ discrepancy 
with the state of the art empirical potential in the re¬ 
gion for which vibrational energies had been measured. 
The ab initio potential also predicted the existence of 
one more vibrational level than the empirical potential, 
in the cases of BeH and BeD. This was all for the ground 
electronic state X (l^S+), so it is of interest to see how 
well the most state-of-the-art ab initio potential for the 
6 e“ Li 2 excited state 6 (l^n„) will be. 

A Fock space MRCC method based on the (2,0) sec¬ 
tor of the Fock space, called FS-CCSD(2,0) (83], was re¬ 
cently implemented and used to calculate potential en¬ 
ergy curves for many excited states of Li 2 [i 2 | with the 
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Tabl© IV. Parameters defining the recommended MLRp^^^ potential 
for the 6(l^nu)-state of all isotopologues of Li 2 . The MLR model is 
defined in Eqs. I19I37I and has damping parameters s — —2 and p — 0.5. 
The appropriate long-range constants are presented in Table [U] The 
fit’s quality was dd — 0.942. 


MLR|f 


De 12 166 cm-i 
re 2.589 848(23) A 


Po 

-0.022 869 54 

Pe 

-5.586 

Pi 

-2.225 706 

Pr 

7.787 

P 2 

-6.019 341 

Ps 

25.55 

Ps 

-10.672 48 

P9 

35.64 

Pa 

-13.746 79 

PlO 

26.5 

P 5 

-12.632 2 

Pll 

8.4 


ANO-RCC basis set [o^. While in principle possible, 
DKH and BOB corrections have not been made in any 
Li 2 ab initio calculations to date. However, fortunately 
we found in Section HO that the addition of S{r) or 
R{r) functions did not significantly improve the fit to 
the data, meaning that Born-Oppenheimer breakdown 
beyond the effects included from Eqs. and [Til do not 
seem to have a big effect in this particular state of Li 2 , 
at least in the data region. Said another way, the ab ini¬ 
tio Born-Oppenheimer potential is expected to give good 
predictions of the energies of ®’®Li 2 and ^’^Li 2 , with mass- 
dependent differences accounted for only by the Hamil¬ 
tonian’s kinetic energy operator, as was the case with the 
empirical MLR potential. 

Using the ab initio Born-Oppenheimer potential pro¬ 
vided to us by the authors of [l^, and the MLR poten¬ 
tial described by Table lIVi we used LEVEL to calculate 
the vibrational energies of both the ®’®Li 2 and ^’^Li 2 iso¬ 
topologues. We found that the highest levels had outer 
classical turning points of several thousand Angstroms, 
and therefore we found it useful to use the recently de¬ 
veloped mapping which allows the radial mesh to extend 
to r = oo when numerically solving the Schroedinger 
equation MM. which is also implemented in LEVEL. 
With this method we were able to find up to u = 91 for 
®’®Li 2 and u = 98 for ^’^Li 2 , however, when we calcu¬ 
lated the scattering wavefunction, the number of nodes 
indicated that the highest bound vibrational levels should 
be u = 92 and u = 100 respectively. Impressively, these 
results were identical whether we used the ab initio po¬ 
tential, or the MLR potential. 

We used Le Roy-Bernstein theory to predict these 
missing levels for each isotopologue: For a potential, 
the powers of the binding energies should be linear 

in V [1^. We used the slope calculated from u = 90 and 
91 for ®’®Li 2 , and the slope calculated from u = 97 and 
98, for predicting the energies of u = 92 and u = 99 lev¬ 
els respectively. We then used the last two points again 
to calculate a new slope for predicting the energies of 


V = 100. Interestingly, both the ab initio potential, and 
the MLR potential predict the existence of a ®’®Li 2 level 
bound by < 8 X 10“® cm“^ (< 3 kHz) and a ^’^Li 2 level 
bound by < 2 X 10“^° cm“^ (< 6 Hz). Using Cg^'/r^ and 
Cg'^'/r^ we get that the outer classical turning points for 
the least bound levels of each isotopologue are predicted 
to be at least 13 000 A and 120 000 A respectively. 

These vibrational energies were then used to calculate 
the zero point energies (ZPEs) and vibrational energy 
spacings which are presented in the table below, along 
with the discrepancy between the ah initio and empiri¬ 
cal potentials. We have compared the vibrational ener¬ 
gies (since these are important for photoassociation ex¬ 
periments) and the vibrational spacings (since these are 
important for experiments involving energy transitions). 
For both ®’®Li 2 and ^’^Li 2 , the discrepancy for the vi¬ 
brational energies is less than 12 cm“^. The agreement 
for the vibrational spacings is much better than for the 
case of BeH discussed in the beginning of this subsection. 
The largest discrepancy for a ®’®Li 2 vibrational spacing 
is < 0 . 8 cm“^ and for ^’^Li 2 cm“^ is < 0 . 6 cm“^. 


III. CONCLUSION 

The motivation for this work was to build a potential 
that could predict high-accuracy vibrational energies for 
®’®Li 2 ( 6 ) in the accessible energy range of the recently 
built high-precision experimental setup which has so far 
been very successful for photoassociation spectroscopy 
of c(l^E+) @ and A (l^S+)-states 0. A similar pho¬ 
toassociation apparatus has recently also been setup by 
Kai Dieckmann’s group to measure energy levels of ultra¬ 
cold ®’®Li 2 electronic states dissociating to the 25 -|- 3P 
asymptote . The best ab initio vs empirical potential 
comparison for Li 2 in the literature [isj] , predicted vibra¬ 
tional levels with a discrepancy of up to 2.04 cm“^ for 
the a(l^E+), which would have simply not been good 
enough for Hnding the levels in this type of experiment. 
The spectroscopic features in this type of experiment are 
typically around 0.000 2 cm“^; and covering 2 cm“^ with 
one-minute measurements and a 0.000 2 cm“^ step size 
would take about 7 days. 

However, the empirical MLR potential of for the c- 
state predicted energies were accurate enough to cut the 
experiment’s duration to under 2 days, since the first level 
in the laser’s range turned out to be predicted correctly to 
within 0.525 cm“^, despite this energy being right in the 
middle of a 5000 cm“^ gap in available experimental data 
to guide the empirical potential. In our table comparing 
the ab initio and empirical MLR energies for the 5-state 
we see that the vibrational energies predicted by the ab 
initio potential are sometimes in > 10 cm“^ disagreement 
with the empirical values in the region where the energies 
have in fact been measured. However, the ab initio seems 
to predict all vibrational energy spacings correctly within 
less than 1 cm“^ which is much better than the result in 
the current best ground state 5e“ BeH study Hug. 
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The reason we are interested in measuring more levels 
of the 6-state with high-precision, is because it is sur¬ 
prising that the best ab initio calculation of the hrst 
Li(2^S') — Li(2^P) interaction term (Ca) is still in vast 
disagreement with the empirically fitted values from the 
studies of the A-state [3, Q and c-state UM: despite 
lithium only having 3e“, and this C 3 value having sig- 
nihcance for atomic clocks Q. Lithium is also expected 
to play a major role in polarizability metrology, since po¬ 
larizability ratios can be measured much more precisely 
than individual polarizabilities and Li is the preferred 
choice for the standard in the denominator of such a ra¬ 
tio [^. But this discrepancy in C 3 limits the accura^ 
of a potential Li-based standard for polarizabilities [J. 
Consolingly, in this study we have found that the 6-state 
is predicted to have levels bound by < 8 x 10“® cm~^ 
(< 3 kHz) which would imply an outer classical turning 
point of > 13 000 A, which is larger than any case in our 
awareness. Since the less bound the level measured, the 
more precisely C 3 can be determined from a fit, these 
extremely weakly bound energies are promising for re¬ 
solving the discrepancy. While the technology to mea¬ 
sure these extremely weakly bound energies may still be 
years away, many of the very high vibrational levels pre¬ 
dicted in our analysis are indeed accessible with today’s 

I 


photoassociation technology. 

The least bound levels for the H-state which have been 
measured have binding energies of £'^=83 ~ 6 cm“^, 
£„=88 ~ 0.7 cm“^ and £^=97 « 0.4 cm“^ for ®'^Li2, 
®’®Li2 and ^’^Li2 respectively, and the least-squares fit to 
the data gave a C 3 value with a 95% confidence limit un¬ 
certainty of about ±8 cm“^/A Q, which is currently 
the most precise experimentally determined oscillator 
strength for any system, by an order of magnitude [^. 
The ab initio and empirical MLR potentials for the 6- 
state compared in this work, both give predictions that 
are in great agreement for energy levels that are several 
orders of magnitude less deeply bound than the least 
deeply bound H-state measurements, making it there¬ 
fore possible to obtain an empirical C 3 value far more 
precise than in [3| • Hopefully, this would resolve the age- 
old discrepancy between experiment and theory for this 
C 3 value, which was hrst measured experimentally by 
Loomis and Nusbaum in 1931 [9^. 

Empirical potentials have recently been built for the 6- 
and H-states of: Rb2 in 2009 a nd again in 2013 9S|, 
NaCs in 2009 KCs in 2010 [100|. RbCs in 2010 [m|, 
Cs2 in 2011 [10211 and NaK in 2015 |l03j |. however, this 
is to our knowledge, the frst empirical potential built for 
the 6-state of Li2. 


Table IV. Comparison of the binding energies, denoted G(rii); zero-point energies (ZPE); and vibrational energy spacings, denoted uji = 
— G(i’i); for ®’®Li 2 and ^’^Li 2 . The last column is the difference between the two columns directly prior. Discrepancies of > 0.5 cm~^ 
are marked by one star (two stars if it was for vibrational level within the data range). Lines with blue font are for unobserved levels, lines with 
bold green font are for unobserved levels which are accessible by the UBC lab.All energies were calculated by the program LEVEL 8.2 with atomic 
masses and NUSE=0, IR2=1. ILR=3, NCN=6, and CNN=G 3 . All numbers were converged with respect to the radial mesh parameters, to at least to the 
number of digits shown. 


2015 Empirical 2014 2015 Empirical 2014 

[Present work] ab initio £^(2015)(2014) [Present work] ab initio ci;(2015) — tj(2014) 




De 

-12166 

-12166 

- 


- 

- 

- 

- 


0 

-11977.92296303 

-11978.33396628 

0.00 


ZPE 

188.07703700 

187.66603 370 

0.41 


1 

-11608.64518318 

-11609.63836960 

0.41 


too 

369.27777985 

368.69559 668 

0.58 


2 

-11244.06766785 

-11245.60583257 

0.99 

** 

UJI 

364.57751533 

364.03253 703 

0.54 


3 

-10884.17443003 

-10886.23056706 

1.54 

** 

UJ2 

359.89323782 

359.37526 551 

0.52 


4 

-10528.95288123 

-10531.55359827 

2.06 

** 

LOs 

355.22154880 

354.67696 879 

0.54 


5 

-10178.39360662 

-10181.59083078 

2.60 

** 

CJ4 

350.55927461 

349.96276 749 

0.60 


6 

-9832.49064407 

-9836.21009569 

3.20 

** 

035 

345.90296255 

345.38073 509 

0.52 


7 

-9491.24176614 

-9495.44752764 

3.72 

** 

016 

341.24887793 

340.76256 805 

0.49 


8 

-9154.64865155 

-9159.33297339 

4.21 

** 

017 

336.59311459 

336.11455 425 

0.48 


9 

-8822.71699888 

-8827.85266438 

4.68 

** 

01 s 

331.93165267 

331.48030 900 

0.45 


10 

-8495.45667267 

-8501.03662308 

5.14 

* 

Ole, 

327.26032621 

326.81604 131 

0.44 


11 

-8172.88194455 

-8178.90484561 

5.58 

* 

Olio 

322.57472812 

322.13177 747 

0.44 


12 

-7855.01184367 

-7861.48242495 

6.02 

* 

0111 

317.87010088 

317.42242 065 

0.45 


13 

-7541.87058935 

-7548.78199962 

6.47 

* 

0112 

313.14125432 

312.70042 533 

0.44 


14 

-7233.48805823 

-7240.83072224 

6.91 

* 

Oils 

308.38253112 

307.95127 738 

0.43 


15 

-6929.90023995 

-6937.67179854 

7.34 

* 

0114 

303.58781828 

303.15892 370 

0.43 


16 

-6631.14965407 

-6639.34414481 

7.77 

* 

0115 

298.75058588 

298.32765 373 

0.42 


17 

-6337.28572765 

-6345.89198438 

8.19 

* 

0116 

293.86392641 

293.45216 042 

0.41 


18 

-6048.36515810 

-6057.36911882 

8.61 

* 

0117 

288.92056956 

288.52286 556 

0.40 


19 

-5764.45230323 

-5773.82634065 

9.00 

* 

Oils 

283.91285486 

283.54277 818 

0.37 


20 

-5485.61964752 

-5495.32152707 

9.37 

* 

Olio 

278.83265571 

278.50481 358 

0.33 


21 

-5211.94838867 

-5221.93208012 

9.70 

* 

0120 

273.67125885 

273.38944 695 

0.28 


22 

-4943.52917649 

-4953.75231611 

9.98 

* 

0121 

268.41921218 

268.17976 401 

0.24 


23 

-4680.46301756 

-4690.88456534 

10.22 

* 

0122 

263.06615893 

262.86775 077 

0.20 


24 

-4422.86233998 

-4433.44560079 

10.42 

* 

0123 

257.60067758 

257.43896 455 

0.16 


25 

-4170.85219392 

-4181.56792748 

10.58 

* 

0124 

252.01014606 

251.87767 331 

0.13 
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26 

-3924.57154815 

-3935.39380033 

10.72 

* 

U32b 

246.28064577 

246.17412 715 

0.11 


27 

-3684.17463083 

-3695.08166990 

10.82 

* 

^26 

240.39691732 

240.31213 043 

0.08 


28 

-3449.83225353 

-3460.80662810 

10.91 

* 

LU27 

234.34237731 

234.27504 179 

0.07 


29 

-3221.73304927 

-3232.76002819 

10.97 

* 

1^28 

228.09920426 

228.04659 991 

0.05 


30 

-3000.08454633 

-3011.14856240 

11.03 

* 

1029 

221.64850293 

221.61146 579 

0.04 


31 

-2785.11398599 

-2796.19367267 

11.06 

* 

1^30 

214.97056035 

214.95488 973 

0.02 


32 

-2577.06877209 

-2588.13740983 

11.08 

* 

1^31 

208.04521390 

208.05626 284 

-0.01 


33 

-2376.21640980 

-2387.24930291 

11.07 

* 

W32 

200.85236229 
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34 
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11.03 

* 
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35 
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10.96 

* 
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36 
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10.83 

* 
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177.67183 548 
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37 
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* 
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* 
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39 
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* 
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-0.39 
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9.68 

* 

^39 
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9.21 

* 
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* 
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8.64 

* 
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* 
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* 
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* 
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7.29 
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* 

45 
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5.78 

* 
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* 

47 
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* 
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* 

48 
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* 
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* 

49 
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3.69 

* 
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* 

50 
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* 

UJ49 
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51 
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2.72 

* 

060 
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52 
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* 
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53 
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2.09 

* 

052 
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-0.21 


54 
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1.88 

* 
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-0.16 


55 
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-0.18 


56 

-119.29925480 
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1.53 

* 
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57 
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* 
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18.70868 473 

-0.16 


58 
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* 
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59 
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60 
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0.96 
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61 
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62 
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-5.21648541 

-5.39528263 
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74 
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75 
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77 
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78 
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79 
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-8.2x10“® 


80 
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81 
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080 

0.08754764 
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-4.6x10“® 


82 
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9.4x10“® 
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-3.3x10“® 


83 

-0.05001277 

-0.05387381 

6 .2x10“® 
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0.03744142 
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-2.3x10“® 


84 

-0.02698408 
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3.9x10“® 


083 

0.02302869 
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-1.6x10“® 


85 

-0.01354795 
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2.3x10“® 
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0.01343614 
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-1.0x10“® 


86 

-0.00620764 
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085 
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87 

-0.00252196 
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086 
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88 

-0.00086788 
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087 
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0.00184486 
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89 

-0.00023358 
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1.4xl0“'‘ 


088 

0.00063430 

0.00072412 

-9.0x10“® 


90 

-0.00004176 
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4.7x10“® 


089 
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0.00022715 

-3.5x10“® 


91 
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92 
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0 
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1 
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0.90 


LOo 
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342.23570605 

0.52 


2 
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3 
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1.87 

** 
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-10 647.11910591 

-10 649.47046247 

2.35 

** 

OJ3 

329.70657292 

330.18379142 

0.48 

-10 320.93332504 

-10 323.83503748 

2.90 


LU4 

325.63542498 

326.18578086 

0.55 

-9 998.73993844 

-10 002.14677658 

3.41 

** 

(^5 

321.68826091 

322.19338661 

0.51 

-9 680.53599058 

-9 684.39645127 

3.86 

** 

LJe 

317.75032530 

318.20394786 

0.45 

-9 366.32133328 

-9 370.63111765 

4.31 

** 

UJy 

313.76533363 

314.21465729 

0.45 

-9 056.09871321 

-9 060.83607551 

4.74 

** 

LOs 

309.79504214 

310.22262007 

0.43 

-8 749.87385050 

-8 755.02510535 

5.15 

** 

CJg 

305.81097016 

306.22486271 

0.41 

-8 447.65555774 

-8 453.21503238 

5.56 

** 

t^io 

301.81007297 

302.21829276 

0.41 

-8 149.45592611 

-8 155.42591892 

5.97 

** 

1^11 

297.78911345 

298.19963163 

0.41 

-7 855.29057880 

-7 861.67464296 

6.38 

** 


293.75127597 

294.16534731 

0.41 

-7 565.17897189 

-7 571.97154020 

6.79 

** 

1^13 

289.70310276 

290.11160691 

0.41 

-7 279.14471358 

-7 286.33906354 

7.19 


W14 

285.63247666 

286.03425831 

0.40 

-6 997.21587514 

-7 004.81028500 

7.59 

** 

1^15 

281.52877854 

281.92883844 

0.40 

-6 719.42527720 

-6 727.41708681 

7.99 

** 

1^16 

277.39319819 

277.79059793 

0.40 

-6 445.81074981 

-6 454.19261209 

8.38 

** 

Wl7 

273.22447472 

273.61452739 

0.39 

-6 176.41537833 

-6 185.17792610 

8.76 

** 


269.01468598 

269.39537148 

0.38 

-5 911.28775836 

-5 920.41199797 

9.12 

** 


264.76592814 

265.12761997 

0.36 

-5 650.48228761 

-5 659.93548890 

9.45 

** 

U>20 

260.47650906 

260.80547075 

0.33 

-5 394.05952269 

-5 403.80140979 

9.74 

** 

U!21 

256.13407911 

256.42276492 

0.29 

-5 142.08662331 

-5 152.07806776 

9.99 

** 

UJ22 

251.72334203 

251.97289938 

0.25 

-4 894.63789792 

-4 904.84167206 

10.20 

** 

1^23 

247.23639569 

247.44872539 

0.21 

-4 651.79545438 

-4 662.17593685 

10.38 

** 

UI24 

242.66573521 

242.84244353 

0.18 

-4 413.64994883 

-4 424.17762178 

10.53 

** 

UJ25 

237.99831508 

238.14550555 

0.15 

-4 180.30141612 

-4 190.95182023 

10.65 

** 

U26 

233.22580155 

233.34853271 

0.12 

-3 951.86015751 

-3 962.61210740 

10.75 

* 

U}27 

228.33971282 

228.44125861 

0.10 

-3 728.44765475 

-3 739.28322072 

10.84 

* 

tt'28 

223.32888668 

223.41250275 

0.08 

-3 510.19747508 

-3 521.10198215 

10.90 

* 

^<^29 

218.18123857 

218.25017968 

0.07 

-3 297.25612707 

-3 308.21675935 

10.96 

* 

<^30 

212.88522280 

212.94134801 

0.06 

-3 089.78382285 

-3 100.78764610 

11.00 

* 

1^31 

207.42911325 

207.47230422 

0.04 

-2 887.95509525 

-2 898.98514722 

11.03 

* 

UJ32 

201.80249888 

201.82872760 

0.03 
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11.03 

* 

1^33 

195.99261984 

195.99588635 

0.00 

-2 502.00028970 

-2 513.01470900 

11.01 

* 

1^34 

189.97781839 

189.95891920 

-0.02 

-2 318.29707659 

-2 329.27183801 

10.97 

* 

1^35 

183.74287099 

183.70321311 

-0.04 

-2 141.08217128 

-2 151.97922220 

10.90 

* 

1^36 

177.29261581 

177.21490531 

-0.08 

-1 970.60062518 

-1 981.36772751 

10.77 

* 

W 37 

170.61149469 

170.48154610 

-0.13 

-1 807.10765720 

-1 817.69937842 

10.59 

* 

1^38 

163.66834908 

163.49296798 

-0.18 

-1 650.86524375 

-1 661.22516252 

10.36 

* 

1^39 

156.47421590 

156.24241345 

-0.23 

-1 502.13726684 

-1 512.20192623 

10.06 

* 

a;40 

149.02323630 

148.72797691 

-0.30 

-1 361.18285956 
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9.71 

* 
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140.95440727 

-0.35 

-1 228.24756934 

-1 237.53719262 

9.29 

* 
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132.93529022 

-0.42 

-1 103.55200378 
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8.79 

* 

i4;43 
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124.69556556 

-0.50 
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8.22 

* 
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-0.57 
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7.59 

* 

1*^45 
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-0.63 
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6.91 

* 
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-0.67 
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* 

LJ47 
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-0.70 
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5.51 

* 
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82.15611337 
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* 

UJ4g 
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-467.44362166 
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* 
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-0.64 
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-0.58 
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* 
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2.73 

* 
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